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In this paper we consider scalar-tensor theories, allowing for both conformal and disformal cou- 
pHngs to a fluid with a generic equation of state. We derive the effective coupling for both background 
cosmology and for perturbations in that fluid. As an application we consider the scalar degree of 
freedom to be coupled to baryons and study the dynamics of the tightly coupled photon-baryon fluid 
in the early universe. We derive an expression for the effective speed of sound, which differs from 
its value in General Relativity. We apply our findings to the p-distortion of the cosmic microwave 
background radiation, which depends on the effective sound-speed of the photon-baryon fluid, and 
show that the predictions differ from General Relativity. Thus, the /^-distortion provides further 
information about gravity in the very early universe well before decoupling. 

PACS numbers: 



I. INTRODUCTION 

In the search for an explanation of the observed accelerated expansion of the Universe, theorists have questioned the 
vahdity of General Relativity. While the cosmological constant is the easiest explanation and in very good agreement 
with observations so far, its magnitude and origin remain a mystery and other explanations are searched for. Theories 
such as /(i?)-theories, massive gravity, galileon models and other interacting dark energy models have been studied 
in much detail and their phenomenological consequences explored. Most of these theories give rise to deviations from 
the ACDM model which can be looked for with future observations and therefore models will be discarded or heavily 
constrained. We refer to [1] and [2] for excellent reviews. 

In this paper we consider scalar-tensor theories of gravity as an extension of General Relativity. These theories have 
been revived recently, either because they are equivalent to some theories mentioned above (for example f{R) models 
can be written in the form of a scalar-tensor theory) or they describe limits of those theories. Here, we consider not 
only conformal couplings of matter to the scalar field but also disformal couplings, as motivated for example from 
galileon models, certain limits of massive gravity and Lorentz-breaking models of gravity [3, 4]. As we will see, the 
effective couplings of matter to the scalar degree of freedom are much more complicated expressions and depend not 
only on the scalar field itself but also on its time variation. Such models motivate time-varying effective couplings. 

Unfortunately, it has recently been shown that disformal couplings are difficult to constrain with local experiments 
(see e.g. [3, 5, 6]). Therefore, as an application, we study the consequences of modifications of gravity in the radiation 
dominated epoch. In particular, we find a generic expression for the sound-speed of the tightly coupled photon- 
baryon fluid in theories with conformal and disformal couplings and calculate the distortion of the cosmic microwave 
background (CMB) caused by the dissipation of acoustic waves. As is well known, the dissipation of acoustic waves 
injects energy into the photons and therefore produces slight deviations from the blackbody spectrum by producing a 
positive chemical potential fj,. The deviations are small; the chemical potential created by the dissipation of acoustic 
waves is of order /i « 10~* in the standard inflationary scenario and is not in violation of the current constraint 
< 9 X 10^^ (for work on these issues, see e.g. [7-11] and references therein). Proposed experiments such as PIXIE, 
however, reach this sensitivity to search for deviations of the order /i « 10^^ [12]. As such, these observations probe 
the primordial power spectrum at very small scales and constrain the inflationary epoch [13, 14]. We point out that 
in general modifications of gravity change the spectral distortions because of the different sound-speed of the coupled 
photon-baryon plasma. For theories in which the field is very heavy the sound-speed is smaller than in General 
Relativity and thus an absence of a fi-type distortion in the CMB spectrum could be explained by a lower sound- 
speed during the epoch when the distortion is created (5 x 10^ < z < 2 x 10^). However, the sound horizon is well 
constrained by measurements of the CMB anisotropics (e.g. the position of the first peak is well known), and, as we 
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will see, this further constrains the type of modified gravity theories which are allowed. It is worth pointing out that 
the /x-distortion is the earliest direct probe of modifications of gravity. 

The paper is organised as follows: in Section 2 we present the action considered in this paper, write down the 
perturbation equations which govern the dynamics of the coupled photon-baryon fluid and derive the effective coupling 
between the scalar field and matter and the effective sound-speed for the coupled photon-baryon fluid. In Section 3 
we calculate the /i— distortion. In Section 4 we present our conclusions and an outlook of future work. 



II. EVOLUTION OF PERTURBATIONS 



A. General Equations 



The action we are considering is of scalar-tensor form, namely 



+ S'mattor(X»,5u2), (1) 



in which TZ is the Ricci scalar, Xi are the matter fields in the theory (relativistic and non-relativistic) , is an additional 
scalar degree of freedom and the metrics g^-^^ are related to the metric g by 

= C^^) ^^)g^^ + ^(«) (0)9^09^0 . (2) 

Each matter species can, in general, couple to a different metric, in which the couplings are described by the functions 
C(0) (called the conformal factor) and D{(f)) (called the disformal factor). However, multiple couplings greatly com- 
plicate the equations of motion and in this paper we will consider a coupling to a single species only. We will return 
to multiple couplings in future work. 

The field equations can easily be obtained from the action above [4]. The scalar field equation is given by 

5^"V^V,0-^+O = O, (3) 



C'.. _ fD , _\ D' 



with 

Q = - V. (^g^.^r^'^ j + ^'/-.A.-^..^'^'' (4) 

and T'^'^ being the energy-momentum tensor of the species coupled to (f), which is consequently not conserved: 

= Q^,^- (5) 

We will write down the equations for the perturbations only and work in the conformal (Newtonian) gauge in which 
the metric is given by 

ds^ = a^ir])[-{l + 2^)dfi^ + 2^)5ijdx'dx^]. (6) 
Perturbing equation ( =) yields the following equations 

k - - (1 + w,)(e, ~ 3$) - 3H f ^ - 5, + ^^5,- ^4 - (7) 

k = -n{i- 3w,)9, - -^9, + k^^ + pI^kH, - fcV, + ^^e, j-^^kHci^, (8) 

l+Wi l + lUi Pi [l+Wi)pi 

whilst the scalar field perturbations obey the Klein-Gordon equation 

5'<j} + iUd^ + {k^ + a^V")5(j) = </>(# + 3$) - 2a^{V' - Qo)* + a^^Q- (9) 
The zero-order part of Q is 

^° " 0/ in^TM 2 1^^^ ' ^^^^ 

2(a^C -I- Dya'^pi — 0^)) 



3 



and the perturbation of Q is 



where 



^ - .r. Z\ I^t^l-^^ + + ^3* + + 65^0], (11) 

a^L + D{a'^pi — 0^) 



62 =3L'0(l + w,), (13) 

63=6i?H0(l + i«O + 2i^02|^^_^^^^^ (14) 

64 = - 3i?H0(l + u;.) - 2i?0 (^^ - ^ + , (15) 

C" fC\^^ CD' D"\ ^ ,2^,1 n'„2„ n'J,2\<9o 



and the subscript i denotes the species the field is coupled to. The expression for Qq agrees with [4] for the case 
Wi = 0. 

For the rest of this paper we shall look at the case when the scalar field is coupled to baryons only. We are treating 
photons and baryons as fluids, coupled via Thomson scattering. For photons and baryons, Eqns. ( ) and ( ) become 



^7= "1^7 + 4*, (17) 

= Ifc^J^ - fc2^^ _|. ;,2^ _^ aneaxidb - &-y), (18) 

-0, + 3a>+^05b-^4--'5Q, (19) 
Pb Pb Pb 



Ob^ - nOb + fc2* + -^-^{9-, - 6b) + —0^6 - —k'5<P , (20) 
R Pb Pb 

where we have added the interaction terms for Thomson scattering and R = Zpb/'^p-f 



B. Tight-Coupling Approximation 

We are interested in scales much smaller than the horizon and on time-scales much smaller than the Hubble 
expansion rate. To derive a second order differential equation for (5-y, we ignore therefore terms which involve the 
Hubble expansion rate, the time-evolution of the background scalar field, and the time-derivatives of the scalar field 
perturbations and the gravitational potential. In this limit, the relevant equations read 



5^^ ~ le^, (21) 

= \k'^5^ - k'^a^ + fc^* + aneariOb - 0^), (22) 

5b^ - Ob, (23) 

Ob = kH + ^^{9, - Ob) ^k'ScP, (24) 
H Pb 

{k'^ + a^V")S(j)=a^SQ, (25) 

where 

^0 = ^TTTTT^^'" 6Q^-—^^—Mb + B,S^l (26) 
2(C + Dpb) a^{C + Dpb) 
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and 



B,=^-Da\V' -Qo), (27) 

B5 = ^ - Dk^ - Da^V" - D'a'V + a^C + i^'p^)^. (28) 
2 Pb 



From equation ('2")) we find 



(29) 



(C + Z?p,)(fc2 + a2-t^")+65Pfc 
To leading order in the tight-coupling approximation r = aneCrT ^ 00 which implies 

0^ « (30) 

and therefore 

^f, « f (31) 

Additionally it can be shown (e.g. [18]) that 



a, = (32) 



and so we can neglect the anisotropic stress. Using equation ( ) on the left-hand side of equation (24) we have 

iie^-eb)^R(^0y-k^-i' + y-k^S(l)^ , (33) 
which can be inserted in equation ( .2) to get 

0^ = Ik^S^ + k^^ + R {e^i - fc^* + ^k^S<f?j . (34) 



Rearranging and using equation ( ) we find 



We can then substitute for Sep using equations (2!)) and (31) and this leads to 

4 

where Cs = l/-\/3(l + R) is the standard sound-speed and 



l^ + {l + 3RT)k^c% = (36) 



{C + Dp,)ik^ + aW')+B5Pb' ^ ' 

From this equation we can read off the modified sound-speed 

di^cUl + iRT), (38) 

which reduces to the expression given in [15] for the purely conformal case. 

To account for Silk damping it is necessary to go beyond leading order in the tight-coupling approximation. Doing 
so (see the Appendix for a derivation) yields 

6^ oc e'^^^e-'''^~^i (39) 

where 



rv 

fs^ Csdrj' 
Jo 



(40) 
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is the sound horizon (which in general differs from its value in General Relativity) and ko is the modified damping 
wavenumber: 



1 _ r°° dz{i + z) 



kl J, 6H{l + R)n,aT 



16 , A o/o , o^T■ , 3(1 + i?) ^2 



15 l + R 



1 - 3(2 + R)T + ' . ' T' 



(41) 



It can be seen that the additional terms are multiplied by i?^/(l + i?). Deep in the radiation dominated epoch when 
Py ^ Pb these terms will be totally insignificant. We found that these modifications are irrelevant for the ^-distortion 
considered in the next section. 

The evolution of Cg depends on the details of the coupling functions C(0) and D((j)) and can be very complicated, 
even in purely conformal theories {D — 0). Potentially, can become negative, signalling instabilities. We will, in 
this paper, not study the consequences of an imaginary sound-speed during the radiation dominated epoch, although 
it would be interesting to study black hole formation on small scales in these scenarios. 

There are several observables which depend on the sound-speed and can therefore be used to search for deviations 
from General Relativity. In this paper we consider two of them, namely the sound horizon at decoupling and the 
/i-distortion. The sound horizon at decoupling determines, for example, the position of the peaks in the anisotropy 
spectrum. Since the position of the first peak is well known, fs{zdcc) can not vary too much from its value in General 
Relativity. As an integral of Cg over time, ^^(zdoc) is dependent mostly on the evolution of Cs for redshifts below 10^ or 
so. The /^.-distortion of the CMB blackbody spectrum is created in the redshift range 5 x 10^ < z < 2 x 10® and probes 
length scales of order fc « 10 Mpc"""^ to A: « 10^ Mpc""'^ . As such, /i does not only provide useful information about 
the primordial curvature perturbation and therefore about about inflationary physics, but also about modifications of 
gravity. In fact, the /i-distortion of the CMB spectrum is the earliest possible direct probe of modifications of gravity 
available to us. We therefore turn our attention to calculate /i in the next section. 



III. /i-TYPE DISTORTION DUE TO DISSIPATION OF ACOUSTIC WAVES 

The evolution of the /^-distortion is given by [7, 8] 

1^-^ + 1.4^, (42) 

where the last term describes the change of /t due to the input of energy into the coupled photon-baryon fiuid and 
the first term describes the thermalization process with ^dc being the double Compton scattering time scale. We will 
give the expression for Q below. 
The solution of this equation is 



/i = 1.4 



t{zi) Pf 



:1.4r..^^e-(^/-)-, (43) 



P7 



1 / Y„ 



2/5 / n h2 \ -2/5 



where 

ZDC = 1.97 X 10® f 1 - - ( ^ ) ) I ) , (44) 

^ V 2 VO.24 J J VO-0224 J ' ^ ' 

and Yp is the primordial helium mass fraction. 

To calculate the energy input, we follow [10, 11] who showed that the energy density of an acoustic wave in the 
photon-baryon plasma can be written as 

2 = P7Y^(^^(-)>P, (45) 

where we here ignore the baryon density (which is much smaller than the energy density in photons p-y), Cs is the 
sound-speed of the wave, 6^ is the photon density contrast, — Pj/p-y ~ | and the average {■■■)p denotes an average 
over one period of oscillation. We have 

/li^k 
(46) 
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where P-yik) is the power spectrum of the photon density contrast. The power spectrum P-y{k) for scales well within 
the horizon can be related to the primordial power spectrum P^{k) by 

P^(fc) = A^(fc)P;(A:), (47) 
where A-y is the transfer function. In General Relativity, the transfer function reads 

A^(fc) = 3cos(fcf,)e-'='/'=°, (48) 

where = / Csdrj is the sound horizon and ku is the diffusion scale. The factor of 3 takes into account for the 
fact that the potential is decaying on super-horizon scales, leading to an enhancement in 5^ [10, 16]. It is here that 
modifications of gravity could play an important role too: if the field is light (to < H) the enhancement might be 
larger or smaller, depending on the coupling functions C and D. On the other hand, in theories such as chameleon 
theories, the mass of the scalar degree of freedom is much larger than the Hubble expansion rate H and the field is 
short ranged. This is the case we are considering in this paper and so we expect the factor of 3 to be a very good 
approximation in the theories we consider here. We will return to the case m < 7? in a separate publication [17] in 
which we discuss the transfer functions for other types of theories of modified gravity. 
The primordial power spectrum can be written [9] 

^2 / i„\ n^-l+^Qln ^ 



P-W»l,45n.L45A,5,(A) (49) 

where — 2.4 x 10^^, fco — 0.002Mpc^^, Ug is the spectral index and a = dus/dlnk is the running of the index. 

In the standard case Cg is independent of the wavenumber fc, but as we have seen, the interaction of the baryons 
with the scalar field causes Cg to be dependent on the wavenumber. We are interested in scales much smaller than the 
interaction range of the scalar field (fc 3> m(a)a) and Cg is only varying very slowly in time. Using equation (45) for 
a wave with a given wave vector k, one can derive the energy density of that wave. Integrating then over all waves 
gives the total energy density 



1 ri^k 

Q = b, I j^dl{k)A^,{k)p;{k) . (50) 



Then using the fact that the average over one oscillation of cos^(x) is 1/2 and treating the photon density as 
effectively constant over the timescales considered we can obtain the energy release per unit redshift 



1.1745 X 10- , , 
Pj J k 



We can then use this expression to evaluate equation (43). 

For the models considered here the effect of modified gravity on /i comes solely from the modified sound-speed. 
How Cg deviates from its evolution in General Relativity depends on the coupling functions C{4>) and D{(j)) and the 
potential V{(l)). Therefore, a plethora of possibilities could be explored. To be specific however, in this paper we focus 
on the purely conformal case for which (with C — e^^"^ and D = 0) the effective sound-speed can be written as 



fc2 + m?o? 

In this case the effective sound-speed is smaller than in General Relativity. Precisely when and how much smaller is 
dictated by how the coupling stength /? and the mass [m? — d'^V/d(p'^) evolve in time. It is instructive to estimate the 
deviation of Cg from its value in General Relativity. We find 

^ ^ m,i3^Rn' ^ 27ni,H^(3^ ^^^^ 

Looking at a regime for which a^m^ ^ k^ and k — 100 Mpc~^, we find A — 10~^ for f3 = 10'^, and A — 0.15 
for (3 = 10^. Therefore, rather large couplings are needed in this case for the sound-speed to deviate significantly 
from its value in General Relativity. As mentioned above, for such a theory to be consistent with the observed CMB 
anisotropics, the sound horizon cannot be modified very much. Therefore A has to become small before decoupling 
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FIG. 1: These graphs show the evolution of the effective sound-speeds, with k = 100 Mpc~^, for two examples. In the left plot 
we have taken b = 7 x 10*, d = 10~^, zq = 10® and rrirat = 350. In the right plot we have taken 6 = 5 x 10"', d = 2 x 10"^, 
zq — 5 X 10^ and rrirat = 350. We also plot the evolution of the standard sound-speed (upper curve in both graphs). 



which requires that either /3 decreases or m increases well before decoupling. In the example that follows we focus on 
the former case for which /3 becomes smaller in time and as a concrete example we consider 



where zq is a redshift after which the coupling becomes rapidly smaller, and b and d are constants. For the mass we 
assume that 



with TOrat being a constant. In Fig. 1 we show the evolution of Cs for a couple of choices for 6, d, Zq and rrirat- In the 
first case, the sound-speed deviates from its value in General Relativity at high redshifts (z > 10^), but approaches 
the standard value quickly for z < 10^. In the second case, the deviation is quite large for z > 400000, so Cg approaches 
Cs at a later time than in the first example and the maximum deviation is greater in this case. In both cases the 
sound-speed approaches its standard value in General Relativity. A motivation for the behaviour of /3 could be that 
General Relativity is an attractor in the radiation dominated epoch (see e.g. [19] for early ideas on such models). 

The corresponding results for the /i-distortion as a function of spectral index Us and running a for the evolution of 
Cs are shown in Fig 2, where we also show the predictions for General Relativity which agree with [13], but with the 
relativistic correction of 3/4 taken into account. As expected, the predicted fi is smaller than in General Relativity 
since Cs is smaller. For the examples shown we see that the predictions for fi are very similar, although the evolution 
of the sound-speed is significantly different. The sound horizon for both cases deviates less than a percent from its 
value in General Relativity. For the examples studied here the effects of the coupling of (j) to baryons on n are of order 
10~^, similar to other, less exotic contributions. Thus the predictions for the /j,-distortion depend on the details of 
the evolution of the coupled photon-baryon fluid and therefore provide a window for modifications of gravity. 

In Fig 3 we show a case for which the deviations of the effective sound-speed from Cs are significant for z > 150000, 
but not below. In this case the predictions for fj, are very similar to in General Relativity, but the prediction for the 
sound horizon deviates by 17 percent. This example is ruled out by observations of the CMB anisotropics and was 
added just for illustration. What becomes clear from these considerations (and from eq. ( !)) is that it is the interplay 
between /3 and m (both their magnitude and evolution) which determines the predictions for /i and fg. Even in the 
case of purely conformal couplings, with different choices for /3 and m a range of possible deviations in either fs or 
fi (or both) can be obtained. If we were to allow for a disformal coupling as well (D ^ 0), the results for ^ and fs 
would also depend on the first derivative of the potential and, in general, on the evolution of 4>, at which point general 
statements about predictions and trends are no longer useful but instead concrete models (i.e. concrete choices for 
C((/)), D{(f)) and V{(f>)) have to be studied [17]. However, the results above show that the /j,-distortion of the CMB is 
a useful tool to constrain modifications of gravity further. 




(54) 



m{z) 



mri,tH{z), 



(55) 
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FIG. 2: The upper graph shows the predictions for the /i-distortion in the case of General Relativity as a function of the 
spectral index ris and the running a. The two lower graphs show the predictions for two examples of modified gravity, with the 
parameters taken as in Fig. 1. Since in our case the sound-speed is smaller than in the standard case, /i is predicted to be 
smaller for a given Us and ct. 

IV. CONCLUSIONS 

The main results of this paper can be summarised as follows: firstly, we have derived the general coupling of a fluid 
(relativistic or non-relativistic) to a scalar field, whose influence is described by the effective metric given in Eq. (2). 
We allowed not only for conformal but also for disformal couplings and the expressions in Section 2 are generic when 
the scalar field is coupled to one species only, whose equation of state is w. As is clear from Section 2, already in 
this case the evolution of the effective coupling (e.g. Qo/ Pi in equation (10) for the background evolution) can be 
rather complicated. The evolution of its perturbation (given in eq. (11)) is even more complicated. Therefore, in those 
scenarios, the coupling is in general a function of time. This opens the door to a rich phenomenology. 

Secondly, focusing on the case of the scalar field coupled to baryons, we have derived the expression for the effective 
sound-speed of the tightly coupled photon-baryon fluid, which differs from the expression in General Relativity. As we 
have pointed out, the /^-distortion of the CMB can be used to constrain the evolution of Cg and therefore constrain 
modifications of gravity at very high redshifts (5 x lO'* < z < 2 x 10^) and small length scales (k > 50 Mpc""'^). 
Therefore, thirdly we have calculated the /i-distortion for a simple case in which the coupling becomes smaller as time 
progresses. In this case, the /i-distortion is smaller than in General Relativity, because the sound-speed is smaller. 
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FIG. 3: In this example we have taken 6 = 2.5 x 10"', d = 10~^, zo = 2.5 x 10^ and m-rat = 350 (and k = 100 Mpc~^ to calculate 
Cs on the left). As can be seen, the sound-speed approaches its standard value only for z < 150000 and the deviation of from 
its value in General Relativity is very small. However, the sound horizon deviates from its value in General Relativity by 17 
percent. This case is not compatible with CMB anisotropics and is shown only for illustration. 

Whether this is generic requires a comprehensive analysis of different choices for the functions C((/)), D{ip) and V{(f>). 
We will study this in future publications [17]. 

The /i-distortion of the CMB spectrum is a useful additional probe for testing gravity and not only for the primordial 
power spectrum of perturbations. As we have seen in this paper, for theories with conformal couplings the coupling 
has to be large (/3 > 10'^) for /i to deviate significantly from its value in General Relativity. A study of disformal 
couplings will be presented elsewhere, in which we also study more generic theories for which the mass of the scalar 
field is small and modifications of the transfer functions could be important [17]. 
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Appendix A: Modifications to Diffusion Damping Scale 

We follow the method outlined in [20]. Damping occurs on very small scales on which the gravitational potentials 
are very small but the photon quadrupole is relevant. We assume that all the perturbation variables vary as e*-^'^'''' 
so that 



6b — and — — (56) 



Equation ('^ ') then becomes 



iujOb ^ -ie^ - 0b) ~ —Ob. (57) 
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Rearranging, we find to second order in f ^ 



1 + 


iuR 








u;2 J 


1 - 


iujR 






T 


LU^ ) 



Inserting this in equation ('^2) gives 



lU! 



3iuj 45f 
and after collecting terms we obtain 



2^2 



^7 

ojR 

T 



1 - 



ujR 

T 



= fc^c^(l + 3i?J^) 



iujk c 



16 3a;2^2 ik'^R^ 



15 P 

We can recognise the first term on the right-hand side as fc^c^, and this allows us to write 



uj = kcg + 
= kcs + 



2tCs 
ik^cl 
2f 



16 Suj^R^ 



15 fc2 

16 i?2 



- 6i?^ J" 



15 1 + i? 



l-3(2 + i?)^+^ii±^^2 



(58) 



(59) 



(60) 



(61) 



where, in the second equality, we have inserted the zero order part of the l/r expansion (wq — kcs) into the first order 
correction. Now using this expression in our ansatz for 5-^ we find 



c ikr^ —k^/k'^, 



with 



dz{l + z) 
6H{1 + R)neCrT 



16 i?2 
15 ^ 1 + i? 



l-3(2 + i?)^+^(i±^^2 



(62) 



(63) 
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